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Abstract 

We present a model for calculating charged-lepton mixing matrices. These matri¬ 
ces are an essential ingredient for predicting lepton flavor-violating rates in the lepton 
number nonuniversal models recently proposed to explain anomalies in R-meson de¬ 
cays. The model is based on work on “constrained flavor breaking” by Appelquist, Bai 
and Piai relating the charged-lepton mass matrix, Aif, to those for the up and down- 
type quarks, M.u,d- We use our recent model of lepton nonuniversality to illustrate the 
magnitudes of flavor-violating B-decay rates that might be expected. Decays with 
final states generally have the highest rates by far. 
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The LHCb Collaboration has reported several features of i?-meson decays involving b —)■ 
transitions that consistently point to a departure from the Standard Model (SM) of 
particle physics: 

• The ratio Rk of the decay rates of —)■ for £ = p, e P 

Rk = ^ = 0-745lo;o?4 (stat) ± 0.036 (syst). (1) 

This result is a 2.6a dehcit from the standard model (SM) prediction, Rk = 1 + 

0(10-^) PElIl]. 

• The direct measurement |5], 

B{B+ ir+/iV")[i,6] = (1-19 ± 0.03 ± 0.06) x 10“^ . (2) 

This is about 30% lower than the SM prediction, B{B^ —)■ K~^= (l.TSlglg) x 

lO-^PEllH]. 

• The observable in B^ —)■ angular distributions exhibits a deficit in two 

bins, quantihed by LHCb as 2.9a for each bin [9]. However, the theoretical error is 

debated [iniiiiiE]. 

These three measurements were made in the low = Mj^ region of 1.0-6.0 GeV^, away from 
charmonium resonances in the spectrum. 

• The joint CMS-LHCb measurement [13] 

B{B, ^ /r+/i-)e.p = (2.81°;6") x = (0.76t°;?°) x B{B, ^ /r+/i-)5M • (3) 

Although this is consistent with the SM prediction [T3|, the central value is about 25% 
low — as it is for Rk [T^ . 

The i?i^-measurement suggests lepton nonuniversality (LNU) occurs in 6 —)■ transi¬ 

tions; the other measurements are consistent in magnitude and sign. It is no wonder, then, 
that they have inspired a number of LNU models of new physics (NP) above the electroweak 
energy scale, involving the exchange of multi-TeV particles [Ml [El HU HU EU EU EU EU EU 

EUEU[23HU[2Ul2U[3U[3Dl32l[33ll3a[35ll36l[37l[3Ul3UllQ|- 

LNU interactions at high energy are accompanied by LFV interactions unless the leptons 
involved are chosen to be mass eigenstates na. Such a choice is an act of fine tuning in the 
absence of a dynamical or symmetry mechanism justifying itj^ Further, since charged leptons 
(and quarks) are massless at Rlnui far above the weak scale, it is difficult to understand the 
motivation or need for flavor-invariant Yukawa couplings there. If the anomalies reported 

^Attempts in this direction are in Refs. [3^ 138) . 
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by LHCb hold up, LFV decays such as B ^ K^*'>fie and B — K^*'> /it should occur at 
rates much larger than in the SM due to tiny neutrino masses alone. The purpose of this 
paper is to present a model for estimating these and other LFV rates implied by new LNU 
interactions. 

LHCb data suggest that LNU is affecting muons but not electrons. To describe this, a 
simple model was adopted in Ref. [U] (hereafter referred to as GGL) in which a heavy Z' 
boson couples only to third-generation quarks and leptons, namely, 

HNP = Gb'^^XfM. (4) 

This chiral structure is consistent with R-decay data which is well £t if the SM and NP 
contributions to the b —)■ s£'^£~ interaction are a product of left-handed currents (LL) [22l 
1231 [211 m]. In Eq. Q, G = = 1/A^p Gp is a new Fermi constant. The primed 

fields refer to the gauge basis, the one in which the charged weak currents are generation- 
universal]^ They are related to mass-eigenstate (unprimed) fields by unitary matrices t/p 
and U[: 

3 3 

h'L = d'L:, = Y.^LdLi, ^ 4.. (5) 

i=l i=l 

The interaction responsible for the discrepancies in Rk, B^ —)■ fjGiji~, Bg —)• and 

Fg is then 


n„p(b ^ = G + H.C.] . (6) 

In GGL we said that the hierarchy of the Gabibbo-Kobayashi-Maskawa (GKM) matrix 
VcKM = for quarks and the apparent preference of the new physics for muons over 

electrons suggest that ^ I^L 32 p — I- Then, in order that this Hamiltonian 

deplete the SM contribution, we assumed GUh^ UX < 0- This sign choice is correct if 
Uf ~ Vckm- In truth, however, we know little about Up and Ff other than Vckm = 
and that plays a similar part in the less well-measured PMNS matrix, Vpmns |13]- 

Experimentalists need better targets. A recent paper by Boucenna, Valle and Vicente [H] 
made a hrst stab at this by assuming Up = Vpmns- In our opinion, Vpmns = U'^Up seems 
likely to be strongly influenced by the unknown neutrino mixing matrix. Ref. [M] discussed 
lepton flavor mixing in the context of LNU due to a leptoquark interaction [22]. While 
we will present results for the Z'-induced Rnpi our scheme for charged-lepton mixing is 
independent of the dynamical nature of LNU and LFV. 

The model we propose for calculating charged-lepton mixing matrices. Up p, is based on a 
recent paper on “constrained flavor breaking” by Appelquist, Bai and Piai (ABP) jUj. They 
assumed a global constrained flavor symmetry group, SU{3)^, broken by just two Yukawa 

^This interaction has been extended in Ref. m to the SU (3) 0 SU (2) 0 C7(l)-invariant form it must have 
ii Affp S> Aew |12|) and used to provide a simultaneous explanation for Rp and A consistent 

gauge model must also be anomaly-free. These extensions are not needed in this paper. 
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Figure 1: Moose diagram for the ABP model of constrained flavor breaking, Ref. [H]. The 
solid links are the input Yukawa matrices chosen by ABP (our case A), the dashed link is 
then predicted. 


spurious. This implies one equation among the three Yukawa matrices, 1^, Y^, Y^. They are 
related to the quark and charged-lepton mass matrices in 

3 

T^mass — ^ ^ Y Y T H.C.^ (7) 

by Aia = Yav/\/2 where v is the Higgs vacuum expectation value. The matrices bring 
these to real, diagonal, positive form: 

Ma = Ma, d^ag = MaU^, {a = U, (!,£). ( 8 ) 


The ABP model is based on the moose diagram of Fig. 1. Requiring (a) that the quark 
doublet and singlet helds, Ql,uji and d^, must be assigned to different 57/(3)’s (to have 
realistic masses and Vckm)] (b) that Ll and likewise be assigned to different 5f/(3)’s; 
and (c) that Ll and en be assigned to 57/(3) groups other than Qls (to avoid oc A4u or 
Aid), leaves six possibilities H. The one chosen by ABP is depicted Fig. Having taken 
Yu and Yi independent. Yd is predicted up to a constant, r]\ 


Yd = r]YuY^ . (9) 

Equivalently, Aid = ^AiuAi\. 

As ABP were interested only in masses of the charged leptons, not their mixing, they did 
not consider assignments differing by an interchange of Ll and cr which just interchanges 55 
and Y^. But, this swaps Aig and Ai\ and, hence, 7/| with 7/^, resulting in different mixing 
factors for LNU/LFV processes. ABP rejected choosing Yu and Yd as the flavor-breaking 
spurious because that implies Yg = rfY^Yu and, hence, unrealistic charged-lepton masses. We 
agree. They also rejected the possibility Y„ = rjYdYg, arguing that it has difficulty obtaining 
a large enough top Yukawa coupling. But r; is a free parameter of unknown origin, and there 
is considerable freedom in choosing the textures for the Yq, so we will consider this case. 
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We will see it is closely related to the case ABP considered. Thus, we consider four cases. 
Written as a relation from which 
//-factor): 

Mf = 


is determined. 

they are 

(ignoring the dimensionful 


(A) 


M-PM, 

(B) 

(10) 

M-/M^ 

(G) 

(M-/ Aljt 

(D) 



By a special choice of quark bases, neither the gauge nor mass bases defined above, 
ABP obtained the mass matrix and, hence, mass ratios for the charged leptons from Eq. 0 
in terms of the known quark masses and CKM matrix elements |13|. No assumption of 
particular quark mass textures was required. While their results are not in agreement with 
data, they are not all that bad, so there is promise in their approach. But, the matrices 
diagonalizing their lepton mass matrix are not transformations from the gauge to the mass 
basis, and thus cannot be used for turning HArp in Eq. Q into predictions of LEV rates. For 
that, we need specific textures for A4u,d in the gauge basis, ones that provide a reasonable 
account of the quark masses and Vckm- 

Fortunately, quark mass textures good enough for our purpose exist; see, e.g.. Ref. [15] . 
We use ones developed in connection with a scenario for solving the strong-CP problem in 
QCD [in] . In this scenario, the phases in M.u,d are rational multiples of vr so that they easily 
satisfy 

6 = arg det Mq = arg det Mu + arg det Md = 0 . (11) 


The mass matrices are Mq=u,d = UlMqW^ where we use MS quark masses renormalized at 
the top-quark pole mass Mt = 173.5 GeV, with eigenvalues Mu = diag(0.00126, 0.611,163.5) GeV 
and Md = diag(0.00264, 0.0522, 2.72) GeVQ Then: 


Mu 


Md 


(0,0) 

(0.01038, -27r/3) 

(0,0) \ 


(0.1325,0) 

(0.5964,0) 

(0,0) 

(12) 

(0,0) 

(0,0) 

(163.5,0) / 


(0.0) 

(0.01112,0) 

(0.01322,0) \ 


(0.01013, IT 

73) (0.05012,0) 

(0.1127,7r/3) 

(13) 

(0,0) 

(0,0) 

(2.721,7r/3) / 



The notation is {\Mq^ij\, arg(Mq^ij)). The motivation for these mass textures is explained in 
Appendix B of Ref. jm). Note that, since quark masses are multiplicatively and universally 
renormalized above M^, the lepton mass textures in Eq. ( [l0| ) are insensitive to QGD running 
from Mt to The GKM matrix obtained by diagonalizing Mu,d, removing its unphysical 

^We also use mi,{mb) = 4.18GeV, rndmc) = 1.275GeV, ms(2GeV) = 95MeV, md{2GeV) = 4.8MeV 
and m„(2GeV) = 2.3MeV [45] . 
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phases and casting it in standard form m , is 


/ (0.976,0) (0.216,0) (0.0045,-0.978) \ 

VcKM = Ul^Ui=\ (0.216, tt) (0.976,0) (0.0415,0) . (14) 

V (0.0075,-0.516) (0.410,3.161) (0.999,0) / 

This reproduces measured CKM matrix entries to within a few per cent, except for Vub and 
Vtd which are within 20% [48] . 

The ^ are obtained (up to a diagonal matrix of pure phases) by diagonalizing 

= . (15) 

For cases A and C (and for cases B and D) 




(16) 


Therefore, the (dimensionless) eigenvalues of are the inverses of those of 

i.e.. 


/ 2 2 2 \ / 
{m^,m mjc = [m 


-2 -2 - 2 \ 

)a: 


(17) 


and Ul j^{C) are the same as Ul ^{A) with their first and third columns interchanged. 
The mass-squared matrices for case A are 




/ (0.1218,0) 
(3.191,-2.50) 
V (3.413,-2.57) 


(3.191,2.50) 

(83.84,0) 

(89.69,-0.0716) 


(3.413,2.57) \ 
(89.69,0.0716) 
(95.98,0) / 


/ (171.3,0) 

(38.40,3.08) 
V (0.1589,1.96) 


(38.40,-3.08) 

(8.617,0) 

(0.03660,-1.11) 


(0.1589,-1.96) \ 
(0.03660,1.11) 
(2.77 X 10-^0) / 


The predicted and measured ratios of the lepton masses are 


melfUr 

= 1.53 X 10■^ 

m^/mr 

= 0.00802 

(cases A, B) 

melmj- 

= 1.53 X 10-^ 


= 0.0191 

(cases C, D) 

me/nrir 

= 2.88 X 10■^ 


= 0.0595 

(Ref. [33]) . 


(18) 


(19) 


( 20 ) 


The predicted ratios are not great, but they do exhibit a qualitatively correct hierarchy. 
Different quark mass textures will lead to different ratios]^ 


"‘We have reproduced the results of Ref. [33] using their Atf. The magnitudes of the corresponding Ufp, 
matrix elements are similar to those in Eqs. (21|22). 
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The lepton mixing matrices for this case are 


Ui{A) = 


UkiA) = 


(0.9808,1.325) 

(0.1515,1.065) 

(0.1231,-2.372) 

(0.02203,2.337) 

(0.1016,2.303) 

(0.9946,-1.325) 


(0.1935,-1.945) (0.02597,3.026) \ 

(0.7149,0.8359) (0.6826,0.5264) , (21) 

(0.6719,-2.369) (0.7304,0.4548)/ 

(0.2176, -0.7759) (0.9758, -0.4548) 

(0.9705,-0.8359) (0.2187,2.627) | , (22) 

(0.1039, -1.321) (0.9062 x lO'^, 1.502) 


where phases have been chosen to make real and positive (see Eq. (15)). The columns 
of Ul{A) are the orthonormal eigenvectors Ve^,v^^,Vr^ of A4^j^j^(A) with rows labeled by 
e', ju', r', and similarly for U^{A). For the hermitian conjugate case, with Aii{B) = M~^A4d, 
the mixing matrix [/[{B) = U^{A). The number of physical phases in depends on the 
nature of the neutrino sector, whether Dirac or Majorana. These phases may induce new 
sources of CP violation in decay, but only by interfering with SM amplitudes. Since LEV 
processes have at most tiny SM amplitudes, their rates involve only absolute values of f/| ^ 
elements. 

These mixing matrices, or ones developed from other quark-mass textures and the ABP 
ansatz, can be used to predict LNU and LEV rates in any NP model of these processes. For 
our Unp, Eq. (|^, the elements of interest in Ul are the third row, Vr'^. In particular, the 
amplitudes for B —)■ and Bg —)■ involve U^*U^y In case A, 32 I ^ 1^^1,331 — 

l/\/2 3> 1^7^311, contrary to our naive expectation that these matrices have a CKM-like 
hierarchy [15]. Even more surprising |D^ 3 i| 3> |f7^32l S> |D^ 33 l- Note that this means that 
the f/| 3 j of case D are CKM-like, as naively expected. These features are a consequence 
of the block-diagonal the latter exhibiting an extreme example of level¬ 

crossing. In turn, these trace back to the textures of A4u and Aid (and the ABP ansatz). 
Aiu is (2 X 2) © (1 X 1) block-diagonal and employs a see-saw to make -C rric without an 
0{mu) matrix element. Approximately the same structure in Aid plus t/]) ~ 1 lead to the 
famous relation tan^c = 6*12 — \/rndJ^s in Vckm- 

Finally, we apply these results to our model “third-generation” Hamiltonian, Eq. (|^, and 
evaluate branching ratios for the LEV processes B —)■ and Bg —)■ ■ Since these 

rates will be proportional to both lepton charge assignments may be combined. 

The hrst order of business is to note that cases B and C are excluded in our model. In those 
cases, 177 ^ 33 / is not much smaller than |77|32/, implying Rk ~ 1 . 

For ^i 7 ^ our model implies (summing over both lepton charge modes) 


B{B 

B{B K(*q+j^-) 

B{B. ^ Iflf) 

B{Bs —)■ p,+/i“) 



2p%p 


mj + 

^ J 


[m, 


m|)^ 


/M 


2m^, 


2p£ 

M, 


Bs 


(23) 

(24) 
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Gase 

R+ ^ K+p^t^ 

B+ K+e^T^ 

B+ K+e^p^ 


A 

D 

Exp. 

1.14 X 10-® 
0.89 X 10-® 

< 4.8 X 10-® 

3.84 X 10“^° 
0.67 X 10-^° 

< 3.0 X 10-5 

0.52 X 10-® 
1.17 X 10-^2 
< 9.1 X 10-® 


Gase 

Bg -)■ p^T^ 

Bg -)• 

Bg -)■ e^p"^ 

Bg — )■ r+r 

A 

D 

Exp. 

1.37 X 10-® 
1.06 X 10-® 

4.57 X 10“^° 
0.80 X 10-^° 

1.73 X 10-^2 
3.91 X 10-^5 
< 1.1 X 10-® 

5.61 X 10“^ 
0.76 X 10-^ 


Table 1: Branching ratios for LFV decays of i?-niesons and Bg —)■ from Eqs. (^22[23 


24), nsing the central valnes of pMP^ of ) ~ (4.29 ± 0.22) x 10 ^ |T] and of 

B{Bs —)■ = (2.8](Q;g) X 10“® [13]. All decays are corrected for phase space (see text). 


Branching ratio limits are from Refs. 


Here [T5] . 


Pnp — 


G Tjd* Tjd 
2 '^L33'^L32 


\U'. 


L32l 


4Gi?T/^*T/ o^EMirrif,) , G jjd* TTd iTTi 
^ -Og + -j 0^33(7^32 I 


= -0.136, 


(25) 


L32I 


TTd* TTd ~ 

Ol33'^L32 — 


where V = Vckm, C'g is the Wilson coefficient for the operator O 9 in 6 —)■ se+e 
V^lVts in the qnark-mass model of Ref. [IHI, and pi is the momentnm of the ontgoing lepton 
in the Bg rest frame. The valne of p^p is obtained from the global-£t resnlt Cg^NP — 
-12% Cg. ,SM m, rather than from Rk alone. This p^p applies to axial-vector amplitndes 
as well because the SM interaction renormalized at mb is pure LL to a good approximation. 

From Eq. (25) and the calculated matrices, one can estimate the G-coupling strength. 
For cases A and D one has G ~ 4.3 x 10“® GeV“^ and 1.8 x 10“® GeV“^. These imply the 
approximate upper bounds Rnp = ^/VG = 4.8 TeV and 745 GeV, respectively. These mass 
scales seem low for a Z', but it must be remembered that it couples primarily to the third 
generation. 

There are two approximations in Eq. (23) as applied to R —>■ ratios. The denom¬ 
inator is best-measured for R+ —)■ it is B{B~^ —)■ = (4.29 ± 0.22) x 10“’^ 

integrated over the full g^-range, 0-22 GeV^ [5]. This integration extrapolates over most of 
the charmonium region, = 8-15 GeV^, ignoring the presence of the narrow resonances. 
Gharmonium resonances do not, of course, influence the numerators, but we do not know 
whether LFV searches will include this region. Second, and potentially more important, LFV 
modes to r’s have smaller phase space than those with p~. For these semileptonic decays, 
this effect is accounted for using the results of Ref. [36]. Employing our own calculations, 
we have corrected for phase space the Bg —>■ decay rates involving r’s. 

Our results are shown in Table As was to be expected for our third-generation Rnp-, 
modes involving pr have the largest rates followed by er with rates smaller by one or more 
orders of magnitude. Rates for the experimentally easier ep modes are very small and may 



















be beyond reach in the near future; the one exception in our model is —)■ in case A. 

The large Bg —)■ pr rate predictions are not yet excluded. The best public limits on these 
LFV modes are also listed in Table [19], 03] . 

In conclusion, it is natural to ask how general are our results; are they to be expected in 
other NP models of the i?-decay anomalies or in other schemes for calculating the mixing 
matrices? Particularly, is the relative importance oi B —)■ X/ir that we found likely to be 
a common feature of such models? It is hard to be sure, of course, but we do believe it is. 
As we emphasized of our "Hatp, the LHCb data strongly points to the third generation, or at 
most just the second and third generations, as the seat of lepton nonuniversality. Further, 
the hierarchy of charged-lepton masses — not unlike that for the quarks — suggests block- 
diagonal mass matrices and, therefore, mixing matrices somewhere along the line from our 
original CKM-like expectation to the ones we found in Eqs. ( 2lp2 ) from the ABP ansatz. 

These expectations can be compared with those obtained within other proposed flavor 
models. Our prediction of a generic enhancement over the SM rate of decay modes involving 
the third generation is also advertised in the class of models discussed in Ref. [3H|, although 
they have unobservable LFV by construction. The only model allowing for a direct compar¬ 
ison is Ref. [35] . For the B ^ K transitions to either er or e/i, the branching-ratio ranges 
predicted in our cases A and D encompass those predicted in their model. In the pr case 
our predictions are above theirs in both A and D cases, although they also predict a relative 
enhancement of this channel with respect to the other LFV modes. 

Therefore, in addition to more hrmly establishing the apparent lepton nonuniversality 
in B decays, it is important that LHCb and other experiments mount searches for lepton 
flavor violation, with special attention to improving signihcantly the limits on pr and even 
er decay modes. 
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